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Context. 

' Aims. We study the clustering properties of Distant Red Galaxies (DRG) to test whether they are the progenitors of local 
, massive galaxies. 

Methods. We use the GOODS-MUSIC sample, a catalog of ~3000 Ks-selected galaxies based on VLT and HST observation of 
the GOODS-South field with extended multi-wavelength coverage (from 0.3 to Sfim) and accurate estimates of the photometric 
redshifts to select 179 DRGs with J — Ks > 1.3 in an area of 135 sq. arcmin. 
O i' Results. We first show that the J — Ks > 1.3 criterion selects a rather heterogeneous sample of galaxies, going from the targeted 
, high-redshift luminous evolved systems, to a significant fraction of lower redshift (1 < 2; < 2) and less luminous dusty starbursts. 
^ ■ These low-redshift DRGs are significantly less clustered than higher-z DRGs. With the aid of extreme and simplified theoretical 
' models of clustering evolution we show that it is unlikely that the two samples are drawn from the same population observed 
, at two different stages of evolution. 
Conclusions. High- 2; DRGs likely represent the progenitors of the more massive and more luminous galaxies in the local Universe 
and might mark the regions that will later evolve into structures of intermediate mass, like groups or small galaxy clusters. 
, Low-2 DRGs, on the other hand, will likely evolve into slightly less massive field galaxies. 
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1. Introduction 

Finding and studying large samples of distant luminous 
and evolved galaxies is fundamental to provide a deeper 
insight on the formation of massive galaxies, a process 
that is commonly perceived as a challenging test for cos- 
mological models of structure formation and evolution. 
For this reason, in the recent past, the study of early type 
galaxies at the highest observable redshifts made use of a 
considerable fraction of large telescope time and occupied 
a substantial part of the astronomical literature. 

The search for passively evolving systems at high 
redshift began with the so-called Extremely Red 
Objects (EROs; see also |Elston, Rieke fc Rieke 1988| 
ICimatti et al. 20021 |McCarthy 2004| ) which reproduce the 
colours of ellipticals at z ~ 1. EROs are relatively "new" 
objects, in the sense that they have been recognized as a 
specific class only around 1990, due to the availability of 
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Near InfraRed (NIR) detectors only at that epoch. Elston, 
Rieke & Rieke (1988) found the first two EROs in a sur- 
vey of 10 sq. arcmin. as resolved galaxies with K ~ 16.5 
and R — K ^ 5. After the optical spectroscopic identifica- 
tions, the two objects turned out to be an evolved galaxy 
at z = 0.8 and a dusty starburst at 2 = 1.44, named 
HRIO. It was clear from this survey and successive in- 
vestigations that the ERO population is heterogeneous in 
its main properties (star formation, mass, age, extinction, 
etc.). 

At present, there are various techniques to find evolved 
galaxies at high-z. ICimatti et al. 19991 utilised the crite- 
rion R — K > 5 (Vega) effective in the redshift interval 
0.8 < z < 1.8 hmited to Ks < 20. |Caputi et al. 2004| 
and lAbraham et al. 20041 used a similar selection / — 
K > 4 (Vega) to select red galaxies with 1 < z < 2. 
IPozzetti fc Manucci 20001 suggested a two colour criterion 
{I — K vs J — K) to separate ellipticals from dusty star- 
bursts at 1 < z < 2, which could be extended at higher 
redshift (2.0 < z < 2.5) using redder colours {J — K 
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vs H — K) . IFranx et al. 20031 proposed a simple pure in- 
frared criterion J — K > 2. 3 (Vega) for z > 2.0. In a 
similar way, ISaracco et al. 20041 selected 3 galaxies with 
J-Ks> 3(Vega) in the HDFS at z > 2.5, plausible can- 
didates for high-z massive galaxies, though the statistic 
is very limited. Recently, IDaddi et al. 20041 suggested to 
isolate early-type galaxies according to the BzK criterion 
[{z -K)ab-{B~ z)ab < -0.2 and {z-K)ab >2.5] effi- 
cient at 1.4 < z < 2.5, and with extension at 2.5 < z < 4.0 
using the RJL colour combination. lYan et al. 20 04 pro- 
posed a new class of objects, the high-z EROS (called 
lEROs) with /^(3.6M)//^(850nm) > 20 to select red 
galaxies at 1.5 < z < 3.0 using MIR data. The physi- 
cal properties of massive galaxies at high-z were also in- 
vestigated bv ISaracco et al. 20051 through spectroscopy of 
a limited sample of massive, evolved galaxies with rela- 
tively bright magnitudes {K < 18.4) at 1.3 < z < 1.7 
on the MUNICS survey. A different approach has been 
adopted for the COMBO- 17 survey, in which the intrin- 
sic colour (U-V) rest frame is utilised to isolate galax- 
ies belonging to the Red Sequence: IBell et al. 20041 used 
the relation {U — V^)rcst > 1-40 — 0.31 • z, efficient at 
0.2 < z < 1.1 according to simulations with spectral 
synthesis code. Finally, [Giallongo et al. 2005| adopted a 
slightly different approach: the bi-modality in {U — V)rest 
is empirically fitted to the observations and could be ex- 
tended up to z ~ 3. 

In this paper we focus on the so-called Distant Red 
Galaxies fDR,Gs: rFranx et al. 2003|l . These galaxies are se- 
lected through a. J — K > 2.3(Vega) criteria, designed to 
be sensitive to galaxies with a large 4000 A break at z > 2. 
IFranx et al. 2003l used this technique in the FIRES survey 
ULabbe et al. 2003|l selecting 14 DRGs in the HDFS, down 
to faint Ks magnitudes {Ks < 24.5 in AB mag). By us- 
ing ultra-deep spectroscopy IvanDokkum et al. 20031 pro- 
vided evidence that the brighter DRGs are indeed galax- 
ies at z ~ 2. Even if the evidence for the existence of old 
and massive galaxies is settled by these observations, how- 
ever the lack of a statistical significant sample of DRGs 
hampered the detailed study of many of their proper- 
ties, in particular their number density, their clustering 
properties and their physical properties like mass, star 
formation, age and spectral energy distribution (SED). 
Recently, |Papovich et al. 2005| have derived a sample of 
153 DRGs from the GOODS South down to a shallower 
limit of Ks — 23.2{AB), with the aim of studying in de- 
tail the specific star formation rate (star formation per 
unit mass star) of DRGs. They found that the bulk of the 
star formation in massive galaxies (M > IO^^Mq) occurs 
at early cosmic epochs and is largely complete by z ~ 1.5. 

Analogously to |Papovich et al. 2005[ we use the ex- 
traordinary dataset provided by the GOODS survey to 
extend these studies. In particular, we will adopt the 
GOODS-MUSIC sample, a i^Ts-selected sample with an 
extended wavelength range (from U to 8.0/Ltm band) that 
we compiled using publicly available data in the Chandra 
Deep Field South region and described at length in 
IGrazian et al. 20061 With this complete sample of DRGs, 



we can define in detail their general properties and refine 
previous investigations bv IFranx et al. 2d()3l which used 
only 14 objects in the FIRES survey, though at a fainter 
magnitude limit. 

The structure of the paper is as follows. In § 2 we 
describe the data used to analyse DRG properties. In § 3 
we select DRGs according to the selection criterion defined 
bv lFranx et al. 20(i3l and provide their number density and 
the redshift distribution. In § 4 we study the clustering 
properties of DRGs selected in the GOODS-South field 
and in § 5 we discuss the link between the DRG population 
at z ~ 2 and the local ellipticals. 

All magnitudes, unless otherwise stated, are given 
in the AB system. A concordance ACDM cosmological 
model = 0.3, = 0.7 and Hq = 70 km s'^Mpc"^) 
has been adopted throughout the paper. 

2. The Data 

We use in this paper the data from the Chandra 
Deep Field South (CDFS; IGiacconi et al. 2000|l . ob- 
tained within the GOODS survey. This is a collabora- 
tion between STScI and ESO l|Renzini et al. 2003|l that 
produced an unprecedented dataset of images, cover- 
ing 135 sq. arcmin. from 0.3 to 8.0/^m down to rel- 
atively faint magnitude limits IjGiavalisco et al. 2004|l . 
In particular, we used the ACS images (release VI. 0, 
IGiavalisco et al. 20041) . the ISAAC database (release VI. 0, 
IVandame et al., in preparation]) and the IRAC datase t 
(release Vl.O enhanced, l&ickinson et al., in preparation| | , 
together with U band photometry from WFI@2.2m ESO- 
MPI and VIMOS reduced by our group. 

Using this public dataset, we have produced a 
high-quality multicolour catalog of galaxies in the 
GOODS-South, that we have named GOODS-MUSIC: 
details about the procedure adopted are discussed in 
IGrazian et al. 20061 We briefly remind here that we have 
used all the publicly available images from U to 8.0 /itm 
{U, B, V, i, z, J, H, Ks, 3.6fi, 4.5^, 5.8fi, 8.0^), in a contigu- 
ous area of 135 sq. arcmin., totalling 14847 objects. In 
particular, to isolate a complete sample of DRGs, we use 
here the iCs-selected sample, that consists of 2931 galax- 
ies. The GOODS survey has a complex, inhomogeneous 
exposure map in the Ks band. To properly derive the sta- 
tistical properties of galaxies in this field, the sample has 
been divided in 6 sub-areas of different magnitude limits, 
as described in details in IGrazian et al. 20061 and in Tab. 
n This information is used in this work when the DRG 
statistical properties are studied, such as their number 
density or clustering properties. The typical magnitude 
limit for most of the sample is about Ks = 23.5, and 
extends down to 23.8 in a limited area. 

In IGrazian et al. 20061 we included spectroscopic in- 
formation for 668 galaxies. Recently, IVanzella et al. 20061 
have released further spectroscopic redshifts in the 
GOODS South region. We used this new release to compile 
a revised sample of 973 galaxies with good spectroscopic 
identification. Out of this number, 815 are in the Ks- 
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Fig. 1. The spectroscopic vs photometric redshifts for 973 
galaxies in the GOODS-MUSIC sample. The accuracy is 
iTz = 0.06 and — 0.03 in the redshift range < z < 6. 

selected sample (~ 28% of the total). For the remaining 
sources, we derived a photometric redshift, as described 
in lGrazian et al. 20061 the redshift accuracy in the range 
< z < 6, as shown in Fig^ on this enlarged spectro- 
scopic sample is cr^ = 0.06, which is the same value pre- 
viously found in lGrazian et al. 20061 If we restrict to the 
340 galaxies with red colours (J — Ks > 0.7), as shown 
m FigEl the redshift accuracy is Cz = 0.08 in the redshift 
range < z < 4. 

Rest-frame physical quantities (such as luminosities, 
mass, age, SFR) are derived by using the synthetic library 
of IBruzual fc Chariot 20031 (hereafter BC03), at the spec- 
troscopic redshift, adopting the same technique already 
described in several previous papers (see Fontana et al. 
2004 for more details). 

3. Selection of DRGs 

3.1. The number density of DRGs 

We have selected DRGs according to the criterion defined 
bv IFranx et al. 20031 (J - Ks > 1.3 in AB system, as 
obtained using the transmission curves for the J and Ks 
filters of ISAAC), which is efficient at z > 2. Fig. O shows 
the effect of this selection criterion of DRGs applied to the 
objects of GOODS-MUSIC sample. 

In the GOODS-South region we find 179 galaxies hav- 
ing J — Ks > 1.3. For the reasons described above, the 
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Fig. 2. The spectroscopic vs photometric redshifts for 340 
red galaxies with J - Ks > 0.7 in the GOODS-MUSIC 
sample. The accuracy is az ~ 0.08 and — 0.05 in the 
redshift range < z < 4. There are only 13 galaxies with 
J — Ks > 1.3 and spectroscopic redshifts (red crosses). 
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Fig. 3. Selection of J - > 1.3 objects in the GOODS- 
South Field. Upper limits in the J band are shown as 
vertical arrows. The horizontal line shows the selection 
criteria for DRGs in the GOODS-South area, while dashed 
line indicates the completeness on the DRG selection due 
to the depth of the J band (26.8 AB at S/N = 1). 
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completeness limit of the survey is not homogeneous, with 
a typical value of Ks — 23.5. We use this sample of DRGs 
to study in particular their number density and their spa- 
tial distribution (clustering). 

The number density of DRGs in the GOODS South 
field is derived through the classical log N — log S distri- 
bution, or the number of objects per sq. arcmin. and per 
magnitude bin in the Ks band. This last quantity is ob- 
tained by following the recipe of lAvni fc Bahcall 1980l 



n{Ks) 



AK 



F7 X! 



i=i 



(1) 



where the sum is on the A^fioid surveys (here, the 
6 areas with different magnitude limits described in 
IGrazian et al. 20061 and in TabQ) and on the A^obj ob- 
jects; Area^^^ represents the accessible area of the j-ih 
survey (this is equivalent to the maximum accessible vol- 
ume when the luminosity function is derived). The DRG 
counts have been computed in bins of AKs = 0.5 magni- 
tude. 

Fig. 01 shows the surface density of DRGs in the 
GOODS South field and compares it with the re- 
sults derived in the HDFS by the FIRES survey 
ULabbe et al. 2003}! . Even if the area of HDFS is smaller 
with respect to the GOODS Survey, the DRG number 
densities in these two independent fields are comparable 
(see also Table Notice, however, that different val- 
ues for the number density of DRGs has been derived 
by using the data of the HDFN l|Lanzetta et al. 19981 
IFontana et al. 20001 IDickinson lOMjl . in which one DRG 
is found at Ks < 23.0, and a limited number at 23 < 
Ks < 24 with upper limit in the J band. The sample 
variance between HDFN and HDFS is due to the limited 
area investigated and stresses the necessity of deriving a 
firm measurement for the number density of DRGs in a 
large and deep survey such as the GOODS-CDFS field. 

3.2. Redshift distribution of DRGs 

The large number of DRGs in the GOODS field makes 
it possible to test the selection criterion and to define 
the window function in redshift for DRGs. Fig. [S] shows 
the distribution of the photometric redshifts of DRGs: the 
spectroscopic sample of DRGs is very limited both in red- 
shifts and Ks magnitudes (only 13 galaxies with 19.7 < 
Ks < 22.9 and 0.65 <z< 3.04). The redshift distribution 
of GOODS DRGs is slightly different from that drawn for 
HDFS bv lFranx et al. 20031 and IDaddi et al. 2003L which 
covers the interval 2 < 2; < 4 with a prominent peak at 
2; ~ 3, and in reasonable agreement with the similar analy- 
sis of Papovich et al 2005. In our GOODS-MUSIC sample 
there are DRGs at lower redshifts (1 < z < 2) with bright 
apparent Ks magnitudes {Ks < 22), which are in prac- 
tice absent in small and deep pencil beam surveys, like 
the HDFS. The redshift distribution clearly shows that 
there is a considerable fraction (77 out of 179, i.e. 43%) 
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Fig. 4. The surface density of selected DRGs in the 
GOODS-South Field (triangles), compared with the esti- 
mate obtained in the HDFS (squares. ILabbe et al. 20031 . 

Table 1. Number density E of DRGs in the Ks band for 
GOODS-South and HDFS fields 



bin 


N 


log(S) 


2 


Ks 


AREA 






mean +la — Icr 






• 2 
arcmin 


20.25 


6 


-1.13 -0.91 -1.38 


1.05 


20.34 


135.372 


20.75 


14 


-0.68 -0.58 -0.82 


1.25 


20.81 


135.372 


21.25 


16 


-0.63 -0.53 -0.75 


1.42 


21.26 


135.372 


21.75 


22 


-0.47 -0.39 -0.57 


1.96 


21.79 


129.692 


22.25 


29 


-0.34 -0.27 -0.43 


2.04 


22.27 


128.273 


22.75 


50 


-0.11 -0.05 -0.17 


2.45 


22.78 


127.935 


23.25 


32 


-0.10 -0.03 -0.19 


2.80 


23.20 


81.272 


23.75 


10 


0.10 0.26 -0.06 


2.75 


23.58 


12.585 


22.50 


4 


-0.13 0.07 -0.30 


2.75 


22.36 


4.500 


23.50 


7 


0.10 0.26 -0.06 


2.75 


23.36 


4.500 


24.50 


8 


0.24 0.38 0.04 


2.75 


24.36 


4.500 


25.50 


18 


0.54 0.65 0.41 


2.75 


25.36 


4.500 



a) the number density E is in units of arcmin^^ mag~^ 

b) bin represents the central bin magnitude in Ks 

c) z and Ks are the mean values of redshift and observed 
Ks magnitude for each magnitude bin 

d) the number density in the second half of the table derives 
from the FIRES survey in the HDFS (|Labbe et al. 2003|l 



of objects at low redshifts {z < 1.9) which satisfy the 
J — Ks selection. With a typical colour J — Ks ~ 1.5, 
they cannot be the result of photometric errors, since this 
should be negligible for relatively bright objects: in fact at 
Ks ~ 21.5 the typical error in magnitude is tr = 0.03. The 
SEDs of these low-redshift DRGs are dominated by power- 
law spectra with a tilt at A ~ 6/zm, which are mostly fitted 
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Fig. 5. Upper panel: the distribution of spectroscopic 
(only 13 objects; short-dashed line) and photometric (solid 
line) redshifts of selected DRGs in the GOODS-South 
Field. The dotted curve is the redshift distribution ob- 
tained for the DRGs using the probability function for 
the redshift for each object derived by the photometric 
redshift code. It is in agreement with the distribution us- 
ing the best estimate for the photometric redshift code. 
The long-dashed line represents the redshift distribution 
for the HDFS f Labbe et al. 2003|l . peaked at z ~ 3. It 
is markedly different from the redshift distribution of the 
GOODS field, since DRGs in the HDFS have fainter Ks 
magnitudes. The redshift distribution of HDFS is compa- 
rable to the redshift distribution of DRGs in the GOODS 
field at Ks > 23 magnitude, as shown in the lower 
panel. Lower panel: the photometric redshift distribution 
for bright {Ks < 22; long-dashed line) and faint {Ks > 22; 
solid line) DRGs. Deep pencil beam surveys (HDFs) pref- 
erentially select objects at z ~ 2, while large area sur- 
veys are biased towards lower- redshift {z < 2) and bright 
{Ks < 22) DRGs (short-dashed line). 

by relatively young galaxies (age/r < 1) and a substan- 
tial amount of extinction {Eb-v ^ 0.5 — 1.0, see Fig. 8 of 
IPapovich et al. 2005| ). 

Fig. El may help in understanding this result, which 
is due to the complex selection effects that are effective 
in this colour criterion. In Fig. we compare the ob- 
served J — Ks colour as a function of redshift with the 
expected J — Ks of a few, selected templates computed 
with the BC03 models. Two of these models are computed 
adopting exponentially declining star-formation histories, 
both started at very high redshift (zform = 20), with solar 
metallicity. The values adopted for the e-folding timescale 
(r = 0.1 and r = 1 Gyr) both produce the same colour at 
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Fig. 6. The J-Ks colour of objects in the GOODS-South 
field as a function of their (spectroscopic or photometric) 
redshift. Upper limits in the J band are displayed as verti- 
cal arrows. The long-dashed horizontal line shows the se- 
lection criteria adopted for DRGs in this paper. The two 
blue solid lines show the J — Ks colour for passively evolv- 
ing galaxies formed at 2; = 20 and with an e-folding star 
formation rate with timescale r = 0.1 and r = 1 Gyr (up- 
per and lower curves, respectively). The red short-dashed 
lines show the same colour for a star-forming galaxy with 
E{B -V) = 1.1 and E{B -V) = 0.5 (upper and lower 
curves, respectively). 

low redshift and show that the J — Ks > 1.3 threshold is 
effective in selecting galaxies at z > 2 that formed their 
stars in a short starburst r < IGyr. At the same time, 
large J — Ks colours may be obtained by star-forming, 
dusty models down to lower redshift z ~ 1. 

This highlights why the DRG population is not a 
unique class of z > 2 objects, but it is contaminated by 
dusty starbursts with z ~ 1.5, whose strong dust absorp- 
tion is responsible for their red infrared colours. The low- 
redshift DRG sub-sample is at the limit of the J — Ks se- 
lection, and can be explained by dust reddening of z ~ 1.5 
star- forming galaxies, as shown in Fig.|Sl If a more drastic 
cut J — Ks colour would be apphed (e.g. J — Ks > 1.8), 
this would ensure a much more efficient selection of galax- 
ies with z > 2, but the sample would be strongly reduced, 
from 179 to 51 galaxies only. 

The difference between low-z and high-z DRG 
has been extensively discussed in a recent paper by 
IPapovich et al. 2005| They argue that lower-z DRG are 
dominated by dusty starbursts, while the higher-z ob- 
jects are made of a more complex stellar population, likely 
a mixture of star-forming, heavily extincted and older, 
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Fig. 7. The distribution of the ratio between the age of 
DRGs and the characteristic timescale r of the exponen- 
tially declining SFR, according to the BC03 spectral syn- 
thesis model. The solid histogram refers to the distribu- 
tion of low-2 DRGs, dominated by relatively young objects 
(age/r < 3) which are typically dusty starbursts, while the 
dashed histogram shows the ratio for 2 < z < 4 DRGs, 
where a considerable fraction (30%) of old and passively 
evolving galaxies arise. 



passively evolving stellar components, with a minority of 
galaxies that are likely genuinely passively evolving. In 
our preliminary, simplified analysis (the most important 
difference with respect to [Papovich et al. 2005| is that we 
do not use models with two-component stellar populations 
and we do not include the 24/i data in the analysis) we 
also have evidence of the same distinction. This is shown 
in Fig[3 where we report the distribution of the ratio be- 
tween the fitted age and the fitted star-formation c-folding 
timescale t (such a ratio is in practice the inverse of the 
Scalo parameter). As it is shown, all low-z DRG are dom- 
inated by actively star-forming, relatively young objects, 
while higher-z DRGs have a broader distribution of age/r, 
including several objects (30% of the high-z DRG sample) 
that are fitted by passively evolving models. 

The average luminosities in the rest-frame I band 
(Vega system) that we infer from the spectral fitting of 
our sample are < M/ -22.3 and < M/ >= -23.2 
at < z >= 1.5 and < z >— 2.7, respectively, and the 
average stellar masses are < >= 8.15 • 10^'^ Mq and 

< >= 9.90 • 10^°Mq (10.76 and 10.88 if we compute 

< log(M) >), respectively. 

It is tempting to speculate on the possible spectral evo- 
lution of these objects. A lower limit to their local lumi- 
nosity can be obtained by assuming that they enter into a 



passive evolution phase soon after we observe them. In this 
case, assuming a truncated star-formation history with so- 
lar metallicity, the BC03 code predicts in the rest-frame I 
band a fading from < z >— 1.5 and < z >— 2.7 to z ~ 
of 2.2 and 2.45 magnitudes, respectively. However, we have 
to take into account that DRG are typically dusty objects, 
such that we should probably normalise this fading to their 
unobscured luminosity. Assuming that the typical redden- 
ing of DRG is E{B - V) ~ 0.75 ± 0.25 with a Calzetti ex- 
tinction curve, and that they evolve to present-day objects 
with little dust extinction, we find that the typical change 
in rest-frame magnitude AM/ = M/(z) — M/(z = 0) 
is 0.26 ± 0.65 at < z >= 1.5 and -0.49 ± 0.65 at 
< z >= 2.7. Given the average rest frame luminosities 
described above, this would imply that the descendents of 
DRG in this simple model have rest frame luminosities 
of about M/(z = 0) = -22.56 and -22.71. The typ- 
ical M* magnitude in the I band for local galaxies in 
the SDSS is Mi = -22.48 (Blanton et al. 2001). which 
increases to —23.2 if one considers only the reddest galax- 
ies {g — r > 0.74). Considering that it is obviously im- 
plausible that all DRGs are observed at the end of their 
star-forming phase, and that therefore they will end up 
in more luminous and massive objects than predicted by 
this exercise, one can conclude that both the low-z and 
high-z DRGs are consistent with being the progenitors of 
local massive galaxies. The analysis of clustering will help 
to clarify this conclusion. 

4. Spatial distribution of DRGs: the clustering 
properties 

It is already known that DRGs are not uniformly dis- 
tributed on the sky, but they are clustered on scales of 
several Mpc. The analysis of the HDFS shows that the 
DRGs are in prevalence concentrated in one quadrant of 
the WFPC, while in the HDFN there is only one DRG, 
suggesting that this population could be strongly clustered 
and affected by cosmic variance l|Vanzella et al. 207)11 
IFranx et al. 20031 IDaHdi et al. 2003jl . such that the small 
area covered by surveys like HDFN or FIRES prevents to 
derive a robust measurement of their clustering properties 
and their redshift evolution. 

We therefore present in the following a detailed anal- 
ysis of the clustering properties of our GOODS-MUSIC 
DRG sample. Thanks to the available statistics, we will 
consider both the overall sample, similarly to what al- 
ready done in previous works, but we will also divide the 
sample into two different sub-groups: the first one, con- 
taining objects with 1 < z < 2, where the dusty starburst 
population is expected to be the dominant component; the 
second one, containing objects with 2 < z < 4, where also 
relatively evolved galaxies are represented in the sample. 

4.1. Angular Two Point Correlation Function 



We have used the 
dLandy fc Szalay 1993| ) to 



Landy-Szalay estimator 
obtain the two-point cor- 
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Fig. 8. The angular distribution of selected DRGs in the 
GOODS-South Field. The symbols are coded according 
to the redshift: DRGs at z > 2 and at z < 2 are shown 
by red triangles and blue circles, respectively; black dots 
refer to normal galaxies at all redshifts. For comparison, 
the size of the HDF is also shown. The DRGs are clustered 
and not uniformly distributed over areas larger than the 
HDFs: this shows that the cosmic variance for DRGs is 
dramatic at small scales. 

relation function (TPCF) in the angular coordinates 
(a,(5) for DRGs in the GOODS Field: 



w{0) = 



GG{9)-2GRi0) + RR{9) 
RR{9) 



(2) 



where GG{9) is the number of observed galaxy pairs at 
distance between 6 and + d9, GR{9) is the number of 
observed-random pairs and RR{9) is the random-random 
pairs. We compute GR and RR as mean values of 1,000 
simulated random catalogs. The random sample of galax- 
ies is obtained by randomly generating the coordinates 
{a,5) in the GOODS-CDFS field. Each random galaxy 
is then retained or rejected according to the magnitude 
limit at the selected position. This ensures to correctly re- 
produce the selection function of observed DRGs, even in 
presence of a complicated exposure map, like the GOODS 
survey one IjGrazian et al. 2006|l . Finally we correct the 
observed w{9) taking into account the bias arising from 
the finite boundary of the sample (see the details in 
Appendix A). 

Errors on the angular correlation function, a^, are de- 
termined by Poisson statistics, through the relation 



'l + w{9) 
GG{9) 



(3) 



We fit the angular correlation function (computed in 
annuli of increasing 9) by a power-law relation, w{9) = 
{9/ 9o)-\ fixing 5 = 0.8. Following 'Gro ft et al. 19971 f see 
also lCroom et al.~ 002 Grazian et al. 2004|l the fit is car- 
ried out by using a Maximum Likelihood Estimator 
(MLE) based on Poisson statistics and unbinned data. A 
detailed description of the MLE can be found in Appendix 
A. 

The results for the DRG TPCF are presented in Fig. 
IHl (large quadrant), together with the MLE fit with the 
corresponding la confidence intervals. Considering the in- 
terval 1 < 9 < 100 arcsec, we find a clustering scale of 
6*0 — 3.19l^;9o arcsec. The mean redshift and absolute 
magnitude for the clustered galaxies are = 2.1 and 
Mj = —22.8, respectively. The small quadrant of Fig. El 
shows the TPCF integrated in circles of increasing aper- 
tures 9. We do not use this quantity to fit the best value 
for since errors are correlated in different bins of an- 
gular separation. However, we can obtain from its value 
an indication of the clustering strength: at = 12 arcsec 
we observe 23 pairs, while simulations of random distribu- 
tions predict 12 pairs only, which is a detection at about 
3cr; at 61 = 6 arcsec we derive an excess of 7 pairs over 3 
random, which represents a 4cr detection. 

By looking at the integrated angular TPCF shown in 
the small quadrant of Fig. |51 we notice that it is still 
significantly non null even at large scales (0 « 50 — 60 
arcsecs), which are comparable to the angular size of the 
HDFs. This result confirms that the difference in the DRG 
number density found in previous surveys is due to both 
the cosmic variance and their strong clustering, whose ef- 
fects can become dramatic when considering deep pen- 
cil beam surveys, which are conducted over small areas, 
like the HDFs or the Hubble Ultra Deep Field (HUDF, 
IReckwith et al. 200311. 

To have a look at the redshift evolution of the DRG 
clustering properties, we compute the correlation scale of 
the low- and high-redshift samples, separately, and find a 
clear evidence of a strong evolution: we indeed estimate 
a correlation scale of 9o = 3.69l3;35 arcsec (z^ ~ 1.5 and 
Mj = —22.30) for the low-redshift sample (76 galaxies), 
and 6*0 = 13.68l^;^g arcsec (z^ = 2.7 M/ = -23.20) for 
the high-redshift one (88 galaxies). 

We note that it is well known that at small scales 
{9 < 10 arcsec) the TPCF is dominated by substructures, 
produced by the existence of multiple galaxies inside mas- 
sive halos (see, e.g. ILee et al. 2005j) . This effect is also 
evident in Fig. |S1 where the presence of close- by galaxy 
pairs or triplets is clearly visible. To measure the clustering 
properties of dark matter halos (DMHs) hosting DRGs, it 
is necessary to avoid using only the smallest scales, where 
the halo occupation distribution (HOD) is plausibly larger 
than unity. Using the total DRG sample, we obtain a cor- 
relation length of 6*0 — 5.89^3;^o arcsec for 6* < 10 arcsec, 
while in the interval 10 < 6* < 100 the TPCF is signif- 
icantly weaker, with a MLE fit of 6*0 = 1.67^150 arcsec 
(see the long dashed lines in Fig. O. It is important to 
remark, however, that the redshift evolution is clearly de- 
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1 2 

log 9(arcsec) 

Fig. 9. Large panel: the differential angular TPCF (in a 
logarithmic scale) for the DRGs in the GOODS field (filled 
circles with Itr error bars). We also plot the MLE best 
fit power-law relation (solid line) with its la confidence 
interval (short-dashed line), as computed in the interval 
1 < 9 < 100 arcsec: the corresponding correlation scale 
of DRGs is 6*0 = 3.19 arcsec. The two dotted lines refer 
to the MLE best fits on limited intervals: 6* < 10 arcsec 
and 6* > 10 arcsec. Note that at small scales the TPCF 
is enhanced by the presence of multiple galaxies in the 
same DMH, while at large scales the boundary effect may 
become critical. Small quadrant: the angular TPCF inte- 
grated over circles of increasing radii (filled circles with Icr 
error bars). 



tected at both scales, although the uncertainties become 
obviously much larger. At the scale oi < 10, indeed, the 
correlation length is 6*0 = 3.84t3'4g arcsec at 1 < z < 2 
and Oq = 15.52t^;^^ arcsec at 2 < z < 4. At 10 < 6* < 100, 
the correlation length is 6*0 = 2.89l;2'g5 ^^rcsec at 1 < z < 2 
and 6*0 = 8.48lg^72° arcsec at 2 < z < 4. Notice that in 
our following discussion we will use the clustering length 
obtained by the fit over the global range 1 < 9 < 100 
arcsec, since it is a robust compromise against boundary 
effects at the largest scales and against HOD effect at the 
smallest scales. 

4.2. Spatial clustering 

To convert the angular correlation length to physical units 
we can invert the angular TPCF through the Limber equa- 
tion IjLimber 1953^ . adopting the DRG redshift distribu- 



tion presented in Fig. Leaving the detailed calculations 
to Appendix A, we have: 

-^obj 

where 1(7) = 3.67909 when 7 = 1.8 is assumed. 

Using the value for 9^) — 3.19l^ gQ derived through the 
MLE fit to the angular TPCF, for the complete DRG sam- 
ple, we obtain a correlation length of tq = 9.781324/1^^ 
Mpc. Using the same Limber equation, the correspond- 
ing comoving correlation lengths are tq = 7.411454/1^^ 
Mpc and rp — 13.36l32oft.~^ Mpc, for the sub-samples at 
1 < z < 2 and 2 < z < 4, respectively. 

We note that the TPCF for the higher-redshift sub- 
sample is different with respect to the value obtained by 
IDaddi et al. 20031 for DRGs in the HDFS, ahhough stiU 
marginally consistent, because of the relatively large error 
budget. We have to notice, however, that for their analy- 
sis they applied a colour selection criterion which is bluer 
{J — Ks > 0.7) than the one adopted here {J — Ks > 1.3). 
We tested that, by selecting in the GOODS region DRGs 
at 2 < z < 4 with their same colour cut, we obtain a sam- 
ple of 232 galaxies, with a typical redshift of z^ = 2.9, hav- 
ing a correlation length of rp = 8.8 ± 1.7h~^ Mpc, which 
is comparable to the value provided bv IDaddi et al. 20031 
(8.3 ± 1.2h'^ Mpc). A redder cut {J -Ks> 1.3) applied 
for DRGs in the HDFS actually results in a larger correla- 
tion length of ro = li.^t^^h-^ Mpc HDaddi et al. 2003|l . 
which is consistent with our estimate. 

As a further comment, we also notice that the error 
associated to our estimate for in our whole sample (^ 
3/i~^ Mpc) is slightly higher than the value quoted by 
IDaddi et al. 20031 for the DRGs in the HDFS, even if the 
samples have a different number of objects (197 DRGs in 
GOODS against 49 in the HDFS). This is due to the fact 
that we include in the error budget the effects of cosmic 
variance, which is the dominant effect in this kind of study 
and it is not included in the error bars quoted for DRGs 
in HDFS. 

It is interesting to compare these results with other 
estimates of clustering strength, for other related classes 
of objects. In order to avoid the dependence of the scale 
length ro on the power-law fit 7, it can be useful to present 
the results in a non-parametric form. This can be done by 
using the quantity ^, defined as the correlation function 
^(r) = {r /ro)^^ integrated over a sphere of a given radius 

^max • 

C(''max) = — / £,{x)x^dx . (5) 

^max 

In general, the larger the scale on which the clustering is 
measured, the easier the comparison with the linear the- 
ory of the structure evolution. Since in the following we 
want to compare our results with those obtained for differ- 
ent values for 7, we prefer to quote clustering amplitudes 
within 20/i~^ Mpc, a scale for which linearity is expected 
to better than a few per cent. Choosing a large radius 
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also reduces the effects of small scale peculiar velocities 
and redshift measurement errors, which can be a function 
of redshift. 

Fig. 1101 compares the values for ^(20ft-~^) that we ob- 
tained for DRGs in the GOODS field (summarised in 
Tabic 12) to the corresponding estimates for other classes 
of objects, both at low and high redshift. It is immedi- 
ately clear that the high-redshift (z > 2) sample of DRG 
is drawn from a remarkably highly clustered population, 
most likely more clustered than the z < 2 DRG popula- 
tion. 

At z = 0, the only galaxies having correlation lengths 
as large as 10-11 h^^ Mpc (corresponding to ^{20h~^) ^ 
1) are morphologically-selected giant ellipticals or radio- 
galaxies. EHii^^^Snil estimate tq — 8.35 ± 0.76/i^^ 
Mpc for early-type galaxies with Mb < —19.5 + 
5log{h) in the Pisces-Perseus super-cluster survey, while 
lAdami &: Mazure 2002l derive a significantly smaller value 
(ro = 7/1-1 Mpc with 7 = -1.79) from the SSRS2 redshift 
survey. The discrepancy between these two measurements 
is probably originated by the presence in the first sur- 
vey of the super-cluster, which enhances the correlation 
function. lOverzier et al. 20031 and [Rottgering et al. 2003| 
find that local radio-galaxies have large clustering lengths 
(see also IPeakcock fc Nicholson 1991|l and that the high 
degree of correlation between hosting ellipticals and lumi- 
nous radio-sources suggests an interesting possible com- 
parison for distant samples. 

For small groups of galaxies in the local 
Universe, the typical value for f has been mea- 
sured by IGirardi et al. 20001 IZandivarez et al. 20031 
IPadilla et al. 20041 and again shown in Fie tTni 
IGollins et al. 200ll report the results of the spatial 
two-point correlation function for the galaxy cluster 
survey _ROSAT-ESO Flux-Limited X-ray (REFLEX), 
finding C(20/i"i) = 2.29 ±0.50 for rich clusters at z < 0.3. 

More ordinary elliptical galaxies show a range of clus- 
tering strength, that is strongly dependent on the abso- 
lute magnitude. We reproduce in Fig. EH the range corre- 
sponding to local elliptical galaxies ranging from Mb = 
-17 to Mb = -2 1, taken from [Norberg et al. 2002| and 
IZehavi et al. 20021 

At high redshifts, we also display the values ob- 
served for EROs IjDaddi et al. 2000|l and for bluer DRG 
IjPaddi et al. 2003|l . 

This compilation of clustering strength for a wide 
range of objects shows that DRG are among the mostly 
clustered objects at galactic scales, and suggests that they 
might be related to the progenitors of similarly clustered 
objects at lower redshifts, as EROs or local massive ellip- 
ticals. Unfortunately, a firm conclusion in this context is 
not straightforward, since we do not know the evolution of 
the bias parameter for this class of high-redshift objects. 
This point will be better discussed in the final section. 



1 1 \ 1 1 1 1 1 \ 1 \ 1 \ I [ 

□ "DRGs (This work) 
OEarly type Galaxies 
^J-Ks>0.7 

^EROs 

X Radio Galaxies 
^ Groups 
^Clusters 




Fig. 10. The integrated clustering strength ^{2Qh~^Mpc) 
as a function of redshifts for different objects: DRGs, 
EROs, powerful radio-galaxies, ellipticals and galaxy 
groups/clusters. Filled squares show the results for low- 
and high-z DRGs in the GOODS region, while void square 
represents the whole DRG sample. The solid lines show 
the predicted evolution of the clustering according to the 
object- conserving model, tuned to the DRGs at low- and 
high-z, while the dashed lines reproduce the clustering 
evolution according to the merging model. The plot sug- 
gests that high-redshift DRGs can be the progenitors of 
local ellipticals, but may evolve into more massive ob- 
jects, like EROs at z ~ 1 and groups/clusters of galax- 
ies in the local universe. The horizontal error bars show 
the redshift intervals for the DRGs in this work. Filled 
triangles show the values of the correlation strength for 
DRGs with J- K s > 0.7{AB) as estimated in the HDFS 
l|Daddi et al. 20 031 at z = 3.1 and in the GOODS region 
at z = 2.9 (this work). 

Table 2. Clustering properties of DRGs 



Type 


ro 


7 




z 


Mi 






Mpc) 










DRGs 


9.78 


+'2.iib 
-3.24 


1.8 


690+"-*"^ 

U.UyU_g ggg 


2.1 


-22.8 


low-z DRGs 


7.41 


+3.45 
-4.84 


1.8 


4iq+0-4i4 

U.^iy_Q 357 


1.5 


-22.3 


high-z DRGs 


13.4 


+2.99 
-3.20 


1.8 




2.7 


-23.2 


J-Ks> 0.7 


8.7711:?^ 


1.8 


0.657_g 272 


2.9 


-23.0 



5. Summary and Discussion 

In this paper we have presented an analysis of Distant Red 
Galaxies (DRG) selected in the GOODS-South region. In 
particular, we have used the GOODS- MUSIC sample, that 
has been compiled from a unique dataset that comprises 
accurate multi-wavelength coverage (14 bands from 0.3 to 
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of ~3000 galaxies in Ks complete sample, with accu- 
rate estimates of the photometric redshifts for all galaxies 
in the field. From the GOODS-MUSIC sample, we have 
selected 179 DRGs according to the criterion proposed 
by IPranx et al. 20031 J — Ks > 1.3 at a typical magni- 
tude Hmit of Ks = 23.5{AB) and down to Ks = 23.8 
in a limited area. The wide and deep covered area (135 
sq. arcmin), together with the extended SED information 
and the precision in photometric redshifts (ctz — 0.06), al- 
lows to study the statistical properties of DRGs, like the 
redshift distribution, number density and clustering prop- 
erties at an unprecedented level. 

The derived number density is consistent with that 
found bv lLabbe et al. 2003L with approximatively 1 DRG 
per sq. arcmin. at Ks = 23.5. The redshift distribution 
shows a smoothed peak around z ~ 2, with extended tails 
both to z = 1 and 2 = 4. Bright DRGs {Ks < 22) 
tend to dominate the 2^1 region, while apparently 
faint DRGs {Ks > 22) are distributed widely around 
z ~ 2.0 — 3.5. The two populations also have different 
intrinsic properties: low-redshift DRG are slightly less lu- 
minous than their higher-z counterparts (< Mj >— —22.3 
and < Ml >— —23.2, respectively), and possibly sHghtly 
less massive (< Mgtar >= 8.15 • lO^^M© and < Mgtar >= 
9.90 • IO^^Mq, respectively). 

In particular, we investigated on the spatial distribu- 
tion of DRGs through the Two-Point Correlation Function 
(TPCF) analysis. We find that DRGs from the overall 
sample are significantly clustered (4(T detection), with a 
typical correlation length of — 3.19t^ gQ arcsec, corre- 
sponding, through the Limber equation and the observed 
redshift distribution, to tq = 9.78tl;llh''^ Mpc. We also 
find that the clustering strength of DRGs increases with 
the J — Ks colour cut used for selection. 

Using the relatively large sample of DRG provided by 
the GOODS-MUSIC sample, we divided the DRG sample 
in two sub-groups in redshift, one with 1 < z < 2 and 
the other with 2 < z < 4. The clustering of low-z DRGs 
is significantly lower than that of the high-z DRGs, with 
ro = 7Altlilh-^ Mpc and UAtlfoh'^ Mpc, respec- 
tively. It is useful to stress here that this behaviour is not 
due to a physical evolution of the DRG population. It is 
the result of a selection criterion which provides an het- 
erogeneous group of dusty star burst and massive/evolved 
galaxies with different redshift distribution. 

Unfortunately, a direct comparison of the clustering 
properties of DRGs with those of other objects can be mis- 
leading, since it is not known a priori the connection be- 
tween these classes. However, it is possible to constrain the 
clustering evolution of the descendents of the DRG popu- 
lation using two extreme, simplified models, as proposed 
by IMatarrese et al. 19971 and IMoscardini et al. 19981 for 
the merging of galaxies in a ACDM hierarchical clustering 
scenario. In one case, that was named object- conserving 
model, we assume that the observed DRGs do not un- 
dergo any subsequent phase of merging with other ob- 
jects, included those of lower mass. This model, which 
is conceptually close to a sort of "passive evolution" sce- 



nario, assumes that the galaxies form at some character- 
istic redshift by some non-linear process which induces a 
bias parameter at that epoch, and that their subsequent 
motion is purely caused by gravity, following the continu- 
ity equation. An obvious consequence of this model is that 
the bias factor will not be constant for all time, but will 
tend to unity as time goes on because the galaxies will be 
dragged around by the surrounding density fluctuations, 
populated by less clustered objects. This scenario, which 
corresponds to have an extremely long merging or disrup- 
tion time, provides an upper limit to the evolution of the 
clustering properties of DRG descendents, and is shown as 
thick solid lines in Fig llOl after normalisations to the DRG 
values obtained in this paper. On the other side, we use a 
merging model, where the - even more extreme - assump- 
tion is that galaxies continue the merging process down 
the lowest redshifts, with the same (high) merger rate of 
their parent halos. This clearly extreme model provides a 
lower limit of the evolution of the clustering properties of 
DRG descendents, and is shown as dashed lines in FietTUI 
These theoretical predictions have been obtained adopt- 
ing the standard ACDM power spectrum, normalised to 
reproduce the local cluster abundance {as = 0.9). 

Although the error budget on the estimate of ^{20h~^) 
on the two DRG samples is still relatively large, we can 
use these two limiting theoretical predictions to attempt 
a physical interpretation of our results. 

First, the observed value of f(20/i~^) for the low-z 
DRGs is outside the range predicted for the evolution of 
the higher z sample: this suggests that is unlikely that 
the two samples are drawn from the same population, ob- 
served at two different stages of evolution. 

If we look at the low redshift range predicted for the 
DRG evolution, it is suggested that high-redshift DRGs 
(i.e. those typically selected at Ks > 22, see FiglS)) likely 
represent the progenitors of the more massive galaxies in 
the local Universe, i.e. the more luminous ellipticals, and 
might mark the regions that will later evolve into struc- 
tures of intermediate mass, like groups or small clusters. 

On the other hand, low-redshift DRGs (i.e. those typ- 
ically selected at Ks < 22), will likely evolve into slightly 
less massive field galaxies, approximately around the char- 
acteristic luminosity L* of local ellipticals. 

Our observations provide further evidence for the so 
called "downsizing" scenario that has emerged in many 
different aspects of high redshift galaxies, providing evi- 
dences that more massive galaxies have formed preferen- 
tially at higher redshifts than less massive ones. Here we 
find the same trend, since high redshift DRGs are more 
clustered, more luminous, and most likely to evolve into 
more massive galaxies than their lower-z counterparts. 
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Appendix A: The Two-Point Correlation Function 
and the Limber Equation 

The calculation of the TPCF over small regions of the sky 
is affected by boundary effects. This bias, known as the in- 
tegral constraint, is produced by the fact that the angular 
TPCF is computed over a limited area fl: the consequence 
is a reduction of the amplitude of the correlation function 

by 



(A.l) 



Following IRoche et al. 20021 we estimated wn numeri- 
cally as 



Y.RR{0)A^e- 



(A.2) 



where we assumed for w{6) a power-law relation: w{9) = 
A^0-^. Fixing S = 0.8 we obtain wn — 10.692^^ and the 
corrected expression for the angular TPCF becomes 



wi0) = w°''%0) + wn 



(A.3) 



The fit of the differential angular TCPF (corrected for 
boundary effects and computed in annuli of increasing 9) 
is carried out by using a Maximum Likelihood Estimator 
(MLE), described in ICroft et al. IMtI This method is 
based on Poisson statistics and unbinned data. Unlike the 
usual minimisation, MLE avoids the uncertainties re- 
lated to the bin size, the position of the bin centre and the 
bin scale (linear or logarithmic). 
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To build this estimator, it is necessary to estimate the 
predicted probability distribution of galaxy pairs, given a 
choice for the correlation length 9o and the slope 5. By 
using all the distances between the random-random pairs 
RR{9), we can compute the number of pairs g{9)d9 in ar- 
bitrarily small bins d9 and use it to predict the expected 
mean number of galaxy-galaxy pairs h{9)d0 in that inter- 
val as 



h{9)d9 = [1 + w{9)]g{9)de , 



(A.4) 



where the correlation function w is modelled by assuming 
a power-law expression, w{9) = {9/9o)~^, S = 0.8. In this 
way, it is possible to use all the distances between the Np 
galaxy-galaxy pairs data to build a likelihood. In partic- 
ular, the likelihood function £ is defined as the product 
of the probabilities of having exactly one pair at each of 
the intervals d9 occupied by the galaxy-galaxy pair data 
and the probability of having no pairs in all remaining 
intervals. Assuming a Poisson distribution, one finds 



c 



Yl exp[~h{9)d9]h{9)d9 J| exp[-h{9)d9] 



(A.5) 



where the index j runs over all the intervals d9 where 
there are no pairs. As usual, it is convenient to define the 
quantity S = — 21n>C, which can be re- written, once we 
retain only the terms explicitly depending on the unique 
model parameter 9o, as 



h{9)d9 - 2J2HH^^)] 



(A.6) 



The integral in the previous equation is computed over 
the range of scales where the fit is made. The minimum 
scale is set by the smallest scale at which we find a DRG 
pair (in our case ^min = 0.6 arcsec), while for the maxi- 
mum scale we adopt 9 max — 15 arcsec. The latter choice 
is made to avoid possible biases from large angular scales, 
where the signal is weak. By minimising S it is possible to 
obtain the best-fitting parameter 6*0. The confidence level 
is defined by computing the increase AS with respect to 
the minimum value of S. In particular, assuming that AS 
is distributed as a with one degree of freedom, AS = 1 
corresponds to 68.3 per cent confidence level. It should be 
noted that by assuming a Poisson distribution the method 
considers all pairs as independent, neglecting their cluster- 
ing. Consequently the resulting error bars can be under- 
estimated (see the discussion in lCroft et al. IQQTjl . 

To convert the TPCF from angular to spatial (3D) co- 
ordinates we can resort to the so-called Limber equation. 
Its original formulation is given by: 



w{9) 



where 



ir '^{n)rldn '^{r2)rU{n2)dr2 



(A.7) 



Adopting the new variables 

ri+r2 n - r2 
r = = (A.9) 

and, assuming the small angle approximation, we obtain, 
ri2=r9{l + yY^\ (A.IO) 
which, when substituted in Ea. llA.7|l . gives 



wi9) 



/,°° ^Hrydr g g[rg(l + yY/^]dy 



(A.ll) 



Finally, using the symmetric properties of ^(r), the ex- 
pression for the Limber equation becomes: 



wi9) = 
where 



0i-''/(7)C*2(r)r^(^)-^dr 



/(7) = 



r(i) 



obj 



3.67909 



(A.12) 



(A.13) 



when the usual value 7 = 1.8 is adopted. 

The redshift distribution of real objects can be written 

as 



dN ^ , , o?r 
= 'firy — 
dz dz 



(A.14) 



where the variation of redshift with comoving distance for 
a AC DM model is given by 



dz _ HoV^Af (1 + zf + »A 
dr c 



(A.15) 



If for the angular correlation function we assume a 
power-law relation w{9) = A9~^ , with (5 = 7— 1, it is easy 
to invert the Limber equation, through the assumption of 
a constant value for rp with redshift: 



,{9) = 



^'-''Iil)rUo°°i^nr{z)V-'<mdz 



(A.16) 



If ^(r) = (7^)"'^ is the TPCF in 3D space, we obtain 
w{9) = (/-)^ for the TPCF in the angular coordinates. 
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